三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率,对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序,该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想,验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合,说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动。为了提高计算规模及计算速度,作者对程序进行了并行化并针对微机机群系统进行了并行优化。实际计算表明本程序具有较高的并行效率。      

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

大气闪烁对近似无衍射激光通信系统的影响

在近似无衍射激光通信系统中,大气湍流的作用导致接收端光强起伏不定,这将对通信系统的误码率产生影响.以光强起伏的概率分布为基础,忽略系统中的其他噪声,推导出折射率结构常数、通信距离和激光波长等影响因素与误码率的关系式.通过计算机模拟,研究了误码率随折射率结构常数的变化规律,讨论了一定湍流强度下通信距离和激光波长对误码率的影响,并以10-9为误码率标准,确定了通信距离的允许范围.结果表明,合理地选择天气和日照条件、减小通信距离以及选择长波长的激光器可以减弱大气湍流对误码率的影响.      

基于相干结构的大涡模拟壁面模型

首先对湍流边界层近壁区相干结构群进行研究,根据共振三波理论使用沿流向周期性发展的多个三维波构造其初始结构,通过直接数值模拟得到该相干结构群的演化特性.在此基础上,提出了基于近壁区相干结构流动参数的新型大涡模拟壁面模型.确定了此壁面模型与主流大涡模拟衔接的原理方案.计算结果表明:该壁面模型在减少计算量的同时,能够为主流区大涡模拟提供重要的壁面大尺度结构信息,使计算取得准确的结果,实现了壁面模型的基本功能.      

运用SIMPLE方法计算大曲率弯道中的湍流流动

在三维任意贴体非正交曲线坐标系下应用同位网格技术,发展了适于求解复杂流动的压力修正算法,采用速度的协变物理分量为动量方程中速度的求解变量,湍流模型选用标准的κ-ε模型.计算表明,方法具有良好的收敛性、准确性和稳定性,计算值与试验数值吻合良好,而且能预示通道内复杂的流动现象.     

固体火箭发动机燃烧室冷态突扩气固多相流动的数值模拟

基于欧拉-欧拉双流体方法,建立固体火箭发动机燃烧室轴对称突扩气固多相流动的统一二阶矩两相湍流模型,采用k-ε双方程模拟气相湍流流动和颗粒相湍流雷诺应力kpg方程描述大尺度颗粒湍流运动、小尺度颗粒各向异性弥散特性以及气固相间作用。数值模拟了突扩流动中的颗粒轴向、径向平均速度和脉动速度以及颗粒-气相脉动速度关联轴向分量的分布特性,计算结果与试验吻合较好。研究结果表明:颗粒轴向速度脉动大于径向速度脉动,轴向平均速度梯度影响脉动速度峰值的发生位置和大小;两相脉动速度关联轴向分量大约是径向分量的2倍。     

用于SiC超光滑表面加工的大气等离子体抛光系统设计

设计了一套基于射频电容耦合放电原理的大气等离子体抛光原型装置,电容耦合等离子体炬首次被应用到了抛光系统中。针对SiC反射镜的加工实验结果表明,系统工作正常,并可以实现亚纳米级的表面粗糙度。     

多段翼型绕流的数值模拟研究

采用C—H型混合的结构化网格,选用Wilcox提出的k-ω模型求解可压的N-S方程,模拟了多段翼型的流场。结果表明：在4°和6°迎角下的非定常计算结果与试验结果的对比,计算结果与试验值吻合的较好。本文应用了亚迭代时间法,正确地描述了30P30N多段翼型前缘结构区域的流动特征和流动机理。     

应用修正的k-ε模型研究超声速H2/Air燃烧

为了研究超声速燃烧中流体可压缩性的影响,对标准k-ε湍流模型进行可压缩性修正(包括结构可压缩性修正和膨胀可压缩性修正两部分)。分别应用标准k-ε模型、修正的k-ε模型和雷诺应力模型(RSM),考虑氢气/空气详细化学反应机理(GR I-M ech 2.11机理,10组分,28基元反应),数值模拟有壁面限制的超声速混合层冷态及热态流场。结果表明:壁面和燃烧对湍流影响都很大;修正模型对冷态以及燃烧场的预测结果优于其它两个;修正模型预测的混合层厚度更薄,燃烧区域更窄,与实验结果吻合地更好。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.